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Abstract 

"Quantum trajectories" are solutions of stochastic differential equations also called 
Belavkin or Stochastic Schrodinger Equations. They describe random phenomena in 
quantum measurement theory. Two types of such equations are usually considered, 
one is driven by a one-dimensional Brownian motion and the other is driven by a 
counting process. In this article, we present a way to obtain more advanced mod- 
els which use jump-diffusion stochastic differential equations. Such models come 
from solutions of martingale problems for infinitesimal generators. These generators 
are obtained from the limit of generators of classical Markov chains which describe 
discrete models of quantum trajectories. Furthermore, stochastic models of jump- 
diffusion equations are physically justified by proving that their solutions can be 
obtained as the limit of the discrete trajectories. 

1 Introduction 

In quantum mechanics, many recent investigations make a heavy use of Quantum Trajec- 
tory Theory with wide applications in quantum optic or in quantum information (cf [14J). 
A quantum trajectory is a solution of a stochastic differential equation which describes 
the random evolution of quantum systems undergoing continuous measurement. These 
equations are called Stochastic Schrodinger Equations or Belavkin Equations (see [7]). 

The result of a measurement in quantum mechanic is inherently random, as is namely 
expressed by the axioms of the theory. The setup is as follows. A quantum system is 
characterized by a Hilbert space TC (with finite or infinite dimension) and an operator 
p, self-adjoint, positive, trace class with Tr[p] = 1. This operator is called a "state" 
or a "density matrix". The measurable quantities (energy, momentum, position...) are 
represented by the self-adjoint operators on H and are called "observable" of the system. 
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The accessible data are the values of the spectrum of the observable. In finite dimension 
for example, if A = Y^=o ^i^i denotes the spectral decomposition of an observable A, the 
observation of an eigenvalue Aj, in the state p, is random and it is obtained with probability: 

P p [to observe AJ = Tr[p P,]. (1) 

Besides, conditionally to the result, the reference state of the system is modified. If 
we have observed the eigenvalue Aj, then the principle called "Wave Packet Reduction" 
imposes the state p to collapse to the new reference state 

i P l pP i . . 

P% TrlpPiY {) 

Quantum Trajectory Theory is then the study of the modification of the state of a system 
undergoing a sequence of measurements. In this way, with the fact that P%P^ — if i j, 
a second measurement of the same observable A, in the state pj, should give 

P p ! [to observe A j] = 1. 

The principle ([2D imposed the new state to be pj = pj. It means that after one measure- 
ment, the information contained in the system is destroyed in the sense that the evolution 
is stopped. 

Actually, in physics applications, a model of indirect measurement is used in order to 
not destroy the dynamic. The physical setup is the one of interaction between a small 
system (atom) and a continuous field (environment). By performing a continuous time 
quantum measurement on the field, after the interaction, we get a partial information of 
the evolution of the small system without destroying it. 

This partial information is governed by stochastic models of Belavkin equations. In the 
literature, there are essentially two different evolutions. 

1. If (p t ) designs the state of the system, then one evolution is described by a diffusive 
equation: 

dpt = L{p t )dt + [p t C* + C Pt - Tr [p t {C + C*)] p t ] dW t , (3) 
where W% describes a one-dimensional Brownian motion. 

2. The other is given by a stochastic differential equation driven by a counting process: 

J(Pt) 



dp t = L(p t )dt + 



Pt 



(dN t -Tr[J(p t )]dt), (4) 



[Tr[J(pt)\ 

where N t is a counting process with intensity j*Tr[J{p s )}ds. 

Equations ([3]) and (TJ1) are called classical Belavkin Equations. The solutions of these 
equations are called "continuous quantum trajectories". Such models describe essentially 
the interaction between a two-level atom and a spin chain f[24].|25j). More complicated 
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models (with high degree of liberty) are given by diffusive evolution with jump described 
by jump-diffusion stochastic differential equations. 

Even in the classical cases (E]) and (jlj), Belavkin equations pose tedious problems in 
terms of physical and mathematical justifications. First rigorous results are due to Davies 
[TO] which has described the evolution of a two-level atom undergoing a continuous mea- 
surement. Heuristic rules can be used to obtain classical Belavkin equations © and (jlj). 
A rigorous way to obtain these stochastic models is to use Quantum Filtering Theory ([7]). 
Such approach needs high analytic technologies as Von Neumann algebra and conditional 
expectation in operator algebra. The physical justification in this way is far from being 
obvious and clear. Furthermore technical difficulties are increased by introducing more 
degrees of liberty and such problems are not really treated. 

A more intuitive approach consists in using a discrete model of interaction called "Quan- 
tum Repeated interactions" . Instead of considering an interaction with a continuous field, 
the environment is represented as an infinite chain of identical and independent quan- 
tum system (with finite degree of liberty). Each part of the environment interacts with 
the small system during a time interval of length h. After each interaction, a quantum 
measurement of an observable of the field is performed. As regards the small system, 
the result of observation is rendered by a random modification of its reference state in 
the same fashion of (EJ). Then the results of measurements can be described by classical 
Markov chains called "discrete quantum trajectories". Discrete quantum trajectories de- 
pend on the time interaction h. By using Markov Chain Approximation Theory (using 
notion of infinitesimal generators for Markov processes), stochastic models for Continual 
Quantum Measurement Theory can be justified as continuous time limit of discrete tra- 
jectories. These models are mathematically justified as follows. Infinitesimal generators 
are obtained as limit (h —>■ 0) of generators of the Markov chains. These limit generators 
give then rise to general problems of martingale Qf9],[I5j). In this article, we show that 
such problems of martingale are solved by solution of particular jump-diffusion stochastic 
differential equations, which should model continuous time measurement theory. This ap- 
proach and these models are next physically justified by proving that the solutions of these 
SDEs can be obtained naturally as a limit (in distribution) of discrete quantum trajectories. 

This article is structured as follows. 

Section 1 is devoted to the description of the discrete model of quantum repeated 
interactions with measurement. A probability space is defined to give account of the 
random character and the Markov chain property of discrete quantum trajectories. Next 
we shall focus on the dependence on h for these Markov chains and we introduce asymptotic 
assumption in order to come into the question of convergence. 

In Section 2, by using Markov chain approximation technics, we obtain continuous time 
stochastic models as limits of discrete quantum trajectories. We compute natural infinites- 
imal generators of Markov chains; these generators also depend on the time interaction h. 
Therefore we obtain infinitesimal generators as limit (h — > 0) of those. It gives then rise to 
general problems of martingale which are solved by jump-diffusion stochastic differential 
equations. 
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Finally in Section 3, we show that discrete quantum trajectories converge in distribution 
to the solution of stochastic differentials equations described in Section 2. The stochastic 
model of jump-diffusion equations is then physically justified as the limit of this concrete 
physical procedure. 

1.1 Discrete Quantum Trajectories 

1.2 Quantum Repeated Measurements 

This section is devoted to make precise the mathematical model of indirect measurement 
and the principle of "Quantum Repeated Interactions". Such model is highly used in 
physical applications in quantum optics or in quantum information (see Haroche [14J). Let 
us start by describing the interaction model without measurement. 

A small system is in contact with an infinite chain of identical and independent quantum 
systems. Each copy of the chain interacts with the small system during a defined time h. 
A single interaction is described as follows. 

The small system is represented by the Hilbert space Ho equipped with the state p. 
A copy of the environment is described by a Hilbert space H with a reference state j3. 
The compound system describing the interaction is given by the tensor product Ho <S> H. 
The evolution during the interaction is given by a self-adjoint operator H tot on the tensor 
product. This operator is called the total Hamiltonian. Its general form is 

H tot = H ®I + I®H + H int 

where the operators Hq and H are the free Hamiltonian of each system. The operator Hi nt 
represents the Hamiltonian of interaction. This defines the unitary-operator 

JJ _ gihHtot 

and the evolution of states of Ho <S> H, in the Schrodinger picture, is given by 

p i-> UpU*. 

After this first interaction, a second copy of H interacts with Ho in the same fashion and 
so on. 

As the chain is supposed to be infinite, the Hilbert space describing the whole sequence 
of interactions is 

T = Ho®<^)H k (5) 

fe>i 

where Hk denotes the /c-th copy of H. The countable tensor product (£) k>1 Hk means 
the following. Consider that H is of finite dimension and that {Xo, X±, . . . , X n } is a fixed 
orthonormal basis of H. The orthogonal projector on CXo is denoted by |Xo)(Xo|. This 
is the ground state (or vacuum state) of H. The tensor product is taken with respect to 
Xq (for details, see [3]). 
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Remark: A vector Y in a Hilbert space 7i is represented by the application \Y) from 
C to 7i which acts with the following way |K)(A) = \XY). The linear form on 7i are 
represented by the operators (Z\ which acts on the vector \Y) by = (Z,Y), where 

( , ) denotes the scalar product of TC. 

The unitary evolution describing the k-th interaction is given by the operator U k which 
acts as U on 7i <S> T~t k , whereas it acts as the identity operator on the other copies of Ti. 
If p is a state on T, the effect of the k-th interaction is: 

p^U kP W k 

Hence the result of the k first interactions is described by the operator V k on B(T) defined 
by the recursive formula: 

Vu+t = U k+1 V k 
Vo — i 

and the evolution of states is then given, in the Schrodinger picture, by: 

P^V kP V k \ (7) 

We present now the indirect measurement principle. The idea is to perform a measurement 
of an observable of the field after each interaction. 

A measurement of an observable of 7i k is modelled as follows. Let A be any observable 
on TC, with spectral decomposition A = J2 P j=o AjPj- We consider its natural ampliation on 

T: ^ 

fc-i 

A k := /. (8) 

j=0 j> k +l 

The result of the measurement of A k is random, the accessible data are its eigenvalues. If 
p denotes the reference state of T, the observation of Xj is obtained with probability 

P[to observe Xj] = Tr[pP 3 k ], j G {0, . . . ,p}, 

where P- is the ampliation of Pj in the same way as (JSj) . If we have observed the eigenvalue 
Xj, the "wave packet reduction" imposes that the state after measurement is 

PfpPf 
93 ~ Tr\pP*} ■ 

Remark: This corresponds to the new reference state depending on the result of the 
observation. Another measurement of the observable A k (with respect to this new state) 
would give P[to observe Xj] = 1 (because P«Pj = if i ^ j). This means that only one 
measurement after each interaction gives a significant information. We recover the phe- 
nomena expressed in the introduction. This justifies the principle of repeated interactions. 

The repeated quantum measurements are the combination of the previous description 
and the successive interactions ((7|). After each interaction, the measurement procedure 
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involves a random modification of the system. It defines namely a sequence of random 
states which is called "discrete quantum trajectory". 
The initial state on V is chosen to be 



i>i 

where p is some state on TLq and each fa = f3 is a fixed state on 7Y. We denote by the 
new state after k interactions, that is: 

Pk = V k p V k *. 

The probability space describing the experience of repeated measurements is fi N * , where 
Q = {0, . . . , p}. The integers i correspond to the indexes of the eigenvalues of A. We endow 
f2 N * with the cylinder cr-algebra generated by the sets: 

Ki,...,i k = W e fi N */^i = h, ■ ■ ■ ,Wfc = 4}- 

The unitary operator C/j commutes with all P fc , for any k and j with k < j. For any set 
{«i, . . . , i k }, we can define the following non normalized state 

/2(ii,...,i fc ) = (/ ® P„ <8> . . . <8> P ifc <8> / . . .) A*jfe (/ (8) Pi! ® . . . (8) Pi fc ® / . . .) 

= (Pi...p>)^(p>...pi). 

It is the non-normalized state which corresponds to the successive observation of the eigen- 
values Aij , . . . , \i k during the k first measurements. The probability to observe these eigen- 
values is 

P[Ai!,...,iJ = P[to observe (X h , . . . , X ik )} = 7V[/2(ii, . . . , i k )\. 

This way we define a probability measure on the cylinder sets of f2 N * which satisfies the 
Kolmogorv Consistency Criterion. Hence it defines a unique probability measure on N . 
The discrete quantum trajectory on T is then given by the following random sequence of 
states: 

~ Pk : ^ — > B(T) 
This next proposition follows from the construction and the remarks above. 



Proposition 1 Let (p k ) be the above random sequence of states. We have for all u G f2 N * : 



Tr 



p^Ut^P^U, 



k+l 



The following theorem is an easy consequence of the previous proposition. 
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Theorem 1 The discrete quantum trajectory (p n ) n is a Markov chain, with values on the 
set of states of Ho <S) {>1 Hi. It is described as follows: 

P [p n+ i = p/ p n = 9 n , . . . , pa = 9 ] = P [p n+ i = p/p n = 9 n \ 

If p n = 9 n , then the random state p n +i takes one of the values: 

4+i 9 n U* +1/ 



P? +1 (Un+^nU^Pr 1 

Tr[(U n+1 6 n U; i+1 )P? +1 ] 
with probability Tr [ (U n+1 9 n U* +1 ) P? +1 ] . 



0,...,p 



In general, one is more interested into the reduced state on the small system H,q only. 
This state is given by taking a partial trace on 7io- Let us recall what partial trace is. If 
7i is any Hilbert space, we denote by Tr^[W] the trace of a trace-class operator W on 7i. 

Definition-Theorem 1 Let Ti, and K be two Hilbert spaces. If a is a state on a tensor 
product 7i <8> K,, then there exists a unique state rj on TC which is characterized by the 
property 

Tr H [r l X}=Tr H ® K [a(X®I)} 

for all X G B(7i). This unique state rj is called the partial trace of a onTi with respect to 
K. 

Let a be a state on T, we denote by E (a) the partial trace of a on TCo with respect 
to <S> k>1 'Hk- We define a random sequence of states on TCo as follows. For all ui in Q N * , 
define the discrete quantum trajectory on Ho 

p n (cu) = E \p n (cu)}. (9) 

An immediate consequence of Theorem 1 is the following result. 

Theorem 2 The quantum trajectory (p n ) n defined by formula (Ej) is a Markov chain with 
values in the set of states on Ho- If p n = Xn, then p n+ i takes one of the values: 



Er 



(J®Pi)l7(xn®/3)l7* {I® P^ 



Tr[U( Xn ®(3)U* (I® Pi)] 
with probability Tr [U (% n ® f3)U* (I ® P<)] . 

Remark: Let us stress that 

(i®Pi)u( X n®[3) u* (i®/*; 



0...p 



Tr[U(xn®P) U* (I® Pi)] 



is a state on Ho ® H. In this situation, the notation E denotes the partial trace on Ho 
with respect to H. The infinite tensor product T is just needed to have a clear description 
of the repeated interactions and the probability space fi N . 
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It is worth noticing that this Markov chain (p^) depends on the time interaction h. By 
putting h = 1/n, we can define for all £ > 



Pnit) = p [nt] . 



(10) 



It defines then a sequence of processes (p n (t)) and we aim to show next that this sequence 
of processes converges in distribution (n — > oo). As announced in the introduction, such 
convergence is obtained from the convergence of Markov generators of Markov chains. The 
following section is then devoted to present these generators for quantum trajectories. 

1.3 Infinitesimal Generators 

In all this section we fix a integer n. Let A be an observable and let p n {t) be the process 
defined from the quantum trajectory describing the successive measurements of A. In this 
section, we investigate the explicit computation of the Markov generator A n of the process 
(p n (t)) (we will make no distinction between the infinitesimal generators of the Markov 
chains (pk) and the process {p n {t)) generated by this Markov chain). For instance, let us 
introduce some notation. 

Let work with T~C = C K+l . The set of operators on 7i can be identified with R p for 
some P (we have P = 2^ +1 ) , we will see later that we do not need to give any particular 
identification). We set E = M p and the set of states becomes then a compact subset of 
IR P (a state is an operator positive with trace 1). We denote by S the set of states and 
E = M p . For any state p G S, we define 



The operators C\ (p) represent transition states of Markov chains described in Theorem 
(12]) and the numbers p l (p) are the associated probabilities. Markov generators for (p n (£)) 
are then expressed as follows. 

Definition 1 Let (p^) be a discrete quantum trajectory obtained from the measurement of 
an observable A of the form A = \P%- Let (p n {t) be the process obtained from (pk) by 
the expression (1 1 . Let define P^ the probability measure which satisfies 




P\P) 



HI Pj) U(n)(p g P)U*{n) (I g Pj) 
°[ Tr[U(n)(p®[3)U*(n)(I®Pi)} 
Tr[U{n){p® P)U*{n)I ® Pi] 



i = 



0...p 



(11) 



P {n) [Pn(0)=p} = l 

P {n) [ Pn {s) = Pk , k/n<s<{k + l)/n] = 1 



(12) 
(13) 
(14) 




where IT n (p, .) is the transition function of the Markov chain (p&) given by 



v 



n n (p,r) = ^p j (p) ( 5 £ ( I1)(p) (r) 



(15) 



i=0 
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for all Borel subset T G B(R P ). 

For all state p G S and all functions f G C%(E) (i.e C 2 with compact support), we 
define 

AJ(p) = nJ(f(p)-f{p))Un{p,dtM) 

= nJ2{f(£ ( r\p))-f(p))p l (p)- (16) 

i=0 

The operator A n is called the "Markov generator" of the Markov chain (or for the 
process (p n (t))). 

The complete description of the generator A n needs the explicit expression of C[ n \p) 
for all p and all i G {0, . . . ,p}. In order to establish this, we need to compute the partial 
trace operation E on the tensor product Tt$®Tt. A judicious choice of basis for the tensor 
product allow to make computations easier. 

Let TCq = C K+1 and let (Qq, . . . , Qk) be any orthonormal basis of TCq. Recall that 
(Xq, . . . , Xtv) denotes an orthonormal basis of TC. For the tensor product we choose the 
basis 

B = (Flo®X ,...,tt K ®Xo,£la®X 1 ,...,tlK®X 1 ,...,tto®X N ,...,tt K ®X N ). 

In this basis, any (N + 1)(K + 1) x (N + l)(K + 1) matrix M on H ® H can be written 
by blocks as a (N + 1) x (N + 1) matrix M = (Mij) Q<i . <N where are operators on Ho- 
Furthermore we have the following result which allows to compute easily the partial trace. 



Claim 1 Let W be a state acting on Ho^H. IfW = (T^ij)o<ij<jv'> the expression ofW 
in the basis B, where the coefficients are operators on Ho, then the partial trace with 
respect to Ti is given by the formula: 

N 

E [W] = J2 W "- 

i=0 

From this result, we can give the expression of the operators C\ n \p). The reference 
state of TC is chosen to be the orthogonal projector on CX Q , that is, with physical notations 

P=\X )(X \. 

This state is called the ground state (or vacuum state) in quantum physics. From general 
result of G.N.S representation in C* algebra, it is worth noticing that it is not a restriction. 
Indeed such representation allows to identify any quantum system (TC, (3) with another 
system of the form (/C, |Xo)(Xo|) where Xq is the first vector of an orthonormal basis of a 
particular Hilbert space /C (see |20j for details). 
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The unitary operator U(n) is described by blocks as U(n) = (Uij(n)) Q<i j <N where the 
coefficients Uy are (K + 1) x (K + 1) matrices acting on TL . For % G {0, . . . , p}, we denote 
Pi = (Pki)o<k,i<N the eigen-projectors of the observable A. Hence the non-normalized states 
Eq[I ® Pi U(h)(p (g) P)U(h)* I ® Pi] and the probabilities p*(p) satisfy 

Eo[/®P^(n)(p®/?)[/(n)*/®Pi] = ^ pL^oHp^oH 

o<k,i<N 

P*(p) = E PL^[^oHp^o(^)]- (17) 

0<k,l<N 

By observing that the operator £ ? - n -'(p) satisfies 

A (p) = ^ , (18) 

for all i G {0, . . . ,p}, we have a complete description of the generator A n . In order to 
consider the limit of A n , we present asymptotic assumption for the coefficient Uij(n) in the 
following section. 

1.4 Asymptotic Assumption 

The choice of asymptotic for U(n) = (t/y(n)) are based on the works of Attal-Pautrat in 
[3J. They have namely shown that the operator process defined for allt > by 

V[ nt ] = U[ nt \(n) ...Ui(n), 

which describes the quantum repeated interactions, weakly converges (in operator theory) 
to a process (Vt) satisfying a Quantum Langevin equation. Moreover, this convergence is 
non-trivian, only if the coefficients Uij(n) obey to certain normalization. When translated 
in our context, it express that there exists operators such that we have for all G 
{0, . . . , N} 2 (recall N + 1 is the dimension of H) 

lim n £ ^ (Uij(n) - SijI) = Ly (19) 

where = |(5oi + <V?)- As the expression (TIT}) given the expression of C[ n \p) only involves 
the first column of U(n), we only keep the following asymptotic 

U 00 (n) = I--L o + o(- 
n \n 

U i0 (h) = -y=L i0 + o f -L ) f or ; > o 



Another fact which will be important in the computation of limit generators is the following 
claim. 
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Claim 2 The unitary condition implies that there exists a self-adjoint operator H such 
that: 



1 - \ 



Furthermore we have for all p G S: 



Tr 



Loop + pLqo + L k opL* ko 



Kk<N 



because Tr[U(n)pU*(ri)] = 1 for all n. 

We can now apply these considerations to give the asymptotic expression of non- 
normalized states and probabilities given by the expression (j!7p . For the non-normalized 
states, we have 



E [J <g> Pi U{n){p ® (3)U{n)* I <g> PJ 
Poo P+ 4^7 XI (Pko L koP + PokP L to) 

v l<fc<JV 

Poo (W + P^oo) + X PliLkopL 



1 

+- 

n 



Kk,KN 



+ ° I - 

n 



(20) 



with probabilities 



p*(p) = Tr[I®P i U{h){p®(3)U{hyi®P l } 



Tr 



Eq[I ® P< [/(ft) (p ® / ® P.-] 



Poo + 



+ -Tr 

n 



X {ploLkop + P l okP L ko) 

<k<N 

Poo ( L ooP + P^oo) + X Pki L kopL 



Kkd<N 



o I - 

n 



The asymptotic expression of C[ n \p) given by the expression (|18|) follows then from 
and fl2T|) . Following the fact that p 00 is equal to zero or not, we consider three cases. 



1. If p 00 = 0, then we have 



4°(p) 



J2l<k,l<NPkl L k0pL^ + o(l) 

Tr[Ei<k,i<NPliLkopL^o + o(l)] ^[£i£W^«Wft]*> 



+ o(l) (22) 
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2. If Pq = 1, then we have 



A (n) (p) 



n 



-TV 



n 



(Loop + pLoo) + Pki L kopLt 

l<k,l<N 

(Loop + pLlo) + Yl Pki L wP L 



Kk.KN 



P + O 



(23) 



3. If Poo ^ {0? 1}; then we have 

A w (p) 



p + ^ 



l 

+- 

n 



rr (Pfeo W + PokP L to) 

"^T Tr f Yl {tiptop + PokpLlo) I x p 
p oo \i<fc<v / 

-T" f Poo (W + P^o) + Pli L koP L m I 

^00 V l<k,l<N J 

+ J^ Tr (^Yl {PwLkoP + pl k pLto)^j xp 

-—TV I p 00 (Loop + P^So) + PkMopL^ ] x p 

Po ° V i<fc,i<JV / 

~Jp^f (Pfco W + PofeP^fco) x Tr f (Pko L koP + PokP L ko) 



Kk<N 



yl<k<N 



+ ° I - 

n 



(24) 



It is worth noticing that all the o are uniform in p because we work in the set of states S 
which is compact. 

With this description, we can now compute the generator limit of A n for any quantum 
trajectory. Next we can establish continuous time model for quantum measurement. This 
is the main subject of the following section. 

2 Jump-Diffusion Models of Quantum Measurement 

In this section, we show that the limit (n — > oo) of generators A n of discrete quantum tra- 
jectories gives rise to explicit infinitesimal generators. From martingale problem technics, 
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we interpret these generators as generators of Markov processes. Besides we show that 
these processes are solution of jump- diffusion stochastic differential equations which are a 
generalization of the classical Belavkin equations (jHJ) and (jlj) presented in Introduction. 

Let us make precise the notion of martingale problem in our framework (see [X9].|15j.[9] 
and [H] for complete references). We still consider the identification of the set of states as a 



compact subset of E = R p for some P. Let II be a transition kernel on E, let a{ 



(•)) 



be a measurable mapping on E with values in the set of positive semi-definite symmetric 
P x P matrices and let b(.) = (&«(.)) be a measurable function from E to E. Let / be any 
C%(E) and let p G E. In this article, we consider infinitesimal generators A of the form 



Af( P ) = 5>(p) 



i=l 



+ 



dpi 



dpidpj 



df(p) ,lf df(p) 

df( P ) 
dpi 



f(p+fi) - f(p) -^2im 



i=i 



n(p, d/i) 



(25) 



The notion of problem of martingale associated with such generators is expressed as follows. 



Definition 2 Let po E E. We say that a measurable stochastic process (pt) on some 
probability space (ft, J 7 , P) is a solution of the martingale problem for (A, po), if for all 

f e cm, 

M{ = f(p t ) - /Oo) - / Af( Ps )ds, t > (26) 
is a martingale with respect to T% = cr(p s , s <t). 

It is worth noticing that we must also define a probability space (Q, J 7 , P) to make explicit 
a solution of a problem of martingale. 

In the following section, we show that Markov generators of discrete quantum trajectory 
converges to infinitesimal generators of the form 025j) . 

2.1 Limit Infinitesimal Generators 

Before to express the proposition which gives the limit infinitesimal generators of A n defined 
in Section 1, we define some functions which appears in the limit. For all % and all state 
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p £ S, set 
9i(p) 



( J2l<kJ<NPkl L kOpL 
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Vi(p) = Tr 
1 

Hp) 
Hp) 



Pki L kopL] 



- P 



J !0 



.Kk.KN 



'Poo 



J2 (Plo L koP + PhkpLto) - Tr [ J2 (Pko L koP + PokpLto)}p 



Kk<N 



Kk<N 



L 00 p + pLqq + L kopLlo- 



(27) 



Kk<N 



This next proposition concerning limit generators follows from results of asymptotic 
described in Section 1. 

Proposition 2 Let A be an observable with spectral decomposition A = ^?=o ^»-^~» where 
Pi = ip\i)o<k,i<N ar & Us eigen-projectors. Up to permutation of eig en-projectors, we can 
suppose that Pq ^ 0. We define the sets 

I = {ie {l,...,p}Mo = 0} and 
J = {l,...,p}\I. 

Let (Pn{t)) be the corresponding quantum trajectory obtained from the measurement of A 
and let A J n be its infinitesimal generator (cf Definition 1). Let A J be the limit generator 
(if it exists) of A J n . It is described as follows. 



1. If I = {1, . . . ,p}, then Poo = 1 and J = ; we have for all f G C^(E): 



lim sup|^/(p)-„4 J /(p)| 

n ^°° pes 



(28) 



where A J satisfies 



A J f(p) = D p f(L(p))+J^ f(p + i^-M-D p f(jM) 



n(p,dp), (29) 



the transition kernel II being defined as 



U(p,dp) = ^2v i (p)S g . {p) {dfj l ). 



i=i 



2. Ifl^{l, 



then Pq 7^ 1 and J ^ ; we have for all f G C^(E): 



lim snp|^/(p)-^l J /(p)| 

n - > °° pes 



(30) 
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where A J satisfies 



A J f(p) = D p f{L{p)) + \ Dlfih^Mp)) 

+ / lf(P + 1*)- f(p) ~ D p f{p)\ n(p, dpi), (31) 
Je 



the transition kernel II being defined as 

U(p,dp) = J2 v i(p) S 9i( P )( d P)- 

Proof: Recall that S is the set of states and it is a compact subset of E. For any 
i G {0, . . . ,p} and for any p G S, let compute 

hm A J J(p) 

n— too 

For this aim, we use asymptotic results of Section 1. As was described, there are three 
cases. 



1. Suppose Poo = 0> we have, 

lim n(f(4 l) (p))-f(p))p\p) 

n— >oo \ / 
f J2l<k,l<NPkl^k0pL^0 

\ T r[Ei<k,i<NPM L kopLt 



Tr 



Yl Pki L kopL* 



.Kk,KN 



(32) 



Moreover, since we have /eC c and since S is compact, the function defined on S 
by 

,f J2i<k,i<N PkiLkopL*o 
\ T r[Y,i<k,i<NPii L kopL^ } 



f(p) 



Tr 



Y Pli L kopL*o 



,Kk,KN 



is uniformly continuous. As a consequence, the asymptotic concerning this case (and 
the fact that all the o are uniform on S cf Section 1) implies the uniform convergence. 

2. Suppose p 00 = 1, by using the Taylor formula of order one, we have 

limn(/(4 t) (p))-/(p))^(p) 



DJ 



(L 00 p + pL* )+ Y Pki L kopL* ] 



l<k,l<N 



-Tr 



(L 00 p + pLlo) + Y Pki L kopLto 



P 



Kk.KN 



(33) 



To obtain the uniform result, we use asymptotic of Section 1 and the uniform conti- 
nuity of Df on S. 
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3. Suppose Pq ^ {0,1}. By applying the Taylor formula of order two, we get the 
convergence 

£ limn(/(£;" ) (p))-/(p))^(p) 

i/P00^{0,l} 



- E 

i/Poo^O- 1 } 



D p f((p 00 ( L ooP + pL* 00 ) + Pki L kopL^o) 

^ l<k,l<N 

-Tr[pi Q (L 00 p + pL* 00 )+ Pli L kopLm]p 

l<k,l<N 

^r D2 p f \ Y {Plo L koP + PokPLto) -Tr[ £ (p l fe0 W + PokP^W ] P 

POO V !< fe <jv l<k<N 

, Y {Pko L koP + PokpLto) - Tr [ Y (Pko L koP + PokpLto) 



Kk<N Kk<N 



(34) 



Let explain more precisely the last equality. When we use the Taylor formula for 
each i such that p l 00 ^ {0, 1}, the term 



P 



Gi(p) = -^=D p fl Y {pio L koP + P ok pLto)-Trl £ {p ko L koP + p okP Ll ) 
appears, but we have 

£ G,(p)=0 

since Ei/ P j *{o,i} Pko = E;/p* £{o,i} Pok = ELoPo* = Y7 i=0 Pko = for any k > 
(indeed we have Ef=o ^» = -0- Furthermore this convergence is uniform for the same 
arguments as previously. 

These three convergence allow us to obtain the two different cases of the proposition. The 
first case of Proposition 2 follows from the first two convergences described above, the 
second case follows from the first and the third convergences above. Before to describe this 
in details, we have to notice that 

p 

Y Y Pki L kopL"lo = Y Lk °P L ko 

i=0 l<k,l<N l<k<N 

since we work with eigen-projectors (Ef=o-^ = -^0- This fact will be used several times. 
Moreover, we have 



Tr[L{p)} = Tr 



Loop + pLq + Y L kopL 



kO 

Kk<N 
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because of Claim 2 in Section 1 concerning the fact that U is a unitary operator. 
Using these facts, in case p 00 = 0, the limit can be written as 



[f(p + P) ~ f(p) ~ Dpf(v)]vi{p)6 gi (j>){diJ.) + D p f{g i {p))v i {p). 



Besides, we have 



D p f(gi(, p))vi{p) = D p f Pki L kopL^o ~ Tr 

\Kk,KN 



Y Pki L kopL*o 

X<k,KN 



Hence it implies the first case of Proposition 2. For / = {1, . . . ,p}, we get indeed 
Aj(p) = D p f(L(p)-Tr[L( P )]p)+ [[f( P + p)-f(p)-D p f(p)]U( P ,dp) 

JE 

= D p f(L(p))+ [ [f( P + p)-f( P )-D p f(p)}Tl( P ,d P ). 

J E 

A similar reasonment gives the expression of the infinitesimal generator in the second 
case where I ^ {1, . . . ,p} and the proposition is proved. □ 

It is worth noticing that generators A 3 are generators of type ( 1251) . it suffices to expand 
the differential terms D p f and D 2 p f in terms of partial derivatives J^- and dp J p - 

In the next section, we present continuous time stochastic models which follows from 
problems of martingale for the limit infinitesimal generators A J . 



2.2 Solutions of Problem of Martingale 

In all this section, we consider an observable A with spectral decomposition 

A = J2^+ Yl V 3 ?' (35) 
iei jeJ(jo 

where I and J are the subsets of {l,...,p} involved in Proposition 2. Let A J be the 
associated limit generator and let po be a state. In order to solve the problem of martingale 
for (A J , po), by Definition[2j we have to define a probability space (Q, J 7 , P) and a stochastic 
process (p J t ) such that the process 

M{ = f(p J t ) - /(po) - f A J f{p J s )ds (36) 

Jo 

is a martingale for the natural filtration of (pi)- 

A classical way to solve the problem of martingale is to define the solution through a 
stochastic differential equation (|H],|12j). 
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Let define a suitable probability space which satisfies the martingale problem. Consider 
(Q, J 7 , P) a probability space which supports a (p + l)-dimensional Brownian motion W = 
(Wo, . . . , W p ) and p independent Poisson point processes (Ni)x<i< p on R 2 and independent 
of the Brownian motion. 

As there are two types of limit generators in Proposition [2, we define two types of 
stochastic differential equations in the following way. Let po be an initial deterministic 
state. 

1. In case J = 0, we define the following stochastic differential equation on P) 

P'l = Po + / L(p J s -)ds + V / / 9i(p J s -)l 0<x<Vl{p J_) [Ni(dx, ds) - dxds] . (37) 

J i=1 JO JR 

2. In case J ^ 0, we define 

Pi = Po + I L(pi_)ds + f hiip^dWiis) 

itJ V{0} ° 

+ / I 9i{p J s -)h <x<Vi ( p J ) [Ni(dx,ds) -dxds]. (38) 

In this way of writing, these stochastic differential equations have a meaning only if 
the process-solution takes values in the set of states (in general the term Vi(p) is not real 
for all operator p). We must modify the expression in order to consider such equation in 
a general way for all process which takes values in operators on 7i - For all i, we define 
when it has a meaning: 

~ , N 12i<k,i<N Ph L kopL*o 

9i{p) = = or f \\ • 

Re(vi(p)) 

Hence we consider the modified stochastic differential equations 

P j t = Po+ L(p J s _)ds + Y^ / / 9i{p J s -) 1 o<x<Re{v,{pi_)) [Ni{dx,ds) - dxds] (39) 
J i=1 Jo Jr 

and 

P J t = Po + I L(p J s _)ds + f K^DdW^s) 

teJUW 

+ Y1 / / ^(^-) 1 o<x<fle(« i ( p J_)) [Ni(dx,ds) - dxds] , (40) 

Let p be a state. The fact that Re(vi(p)) = Vi(p) and g%(p) = g%{p) implies that a solution 
(pi) of the equation ( 1391) (resp (HOI ) is a solution of the equation ( 1371) (resp ( 1381) ) when the 
process (p/) takes values in the set of states. 
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We proceed in the following way to solve the problem of martingale (I36p . Firstly we 
show that the modified equations (l3"9"j) and (j4"Ul) admit a unique solution (we will see below 
that it needs another modification), secondly we show that solutions of ( J39l) and (flOl) 
can be obtained as limit (in distribution) of discrete quantum trajectories (cf Section 3). 
Finally we show that the property of being a process valued in the set of states follows from 
convergence (cf Section 3) and we conclude that solutions of ( |3"9"]) and (BUI) takes values in 
the set of states. Moreover, we show that they are solutions of problem of martingale flHol) . 

The fact that if solutions of (l3"9"l) and fj4*Ul) takes values in the set of states, they are 
solutions of martingale problem ( 1361) is expressed in the following proposition. 

Proposition 3 Let po be any initial state. 

If the modified stochastic differential equations admits a solution (p/) which takes 
values in the set of states, then it is a solution of the problem of martingale (A J , po) in the 
case I = {1, . . . ,p}. 

If the modified stochastic differential equations l \4(ty admits a solution (p[) which takes 
values in the set of states, then it is a solution of the problem of martingale (A J , po) in the 
case J 7^ 0. 

As a consequence, if A J designs the infinitesimal generator of a solution of or 
OPI) . then we have A J f(p) = A J f(p) for all state p and all functions f G C%. 



Proof: Recall we assume that processes take values in the set of states. For any state p, 
we have Re{vi(p)) = Vi(p) and (ji(p) = gi{p) and the part concerning the generators follows. 
Concerning the martingale problem, it is a consequence of the Ito formula. Let p[ = 
(p((t), . . . , Pp(t)) denote the coordinates of a solution of ( 1371) or ( |38l) (with identification 
between the set of operators on Tio and M p ), we have for all /eC c 2 

f(p J t)-f(Po) = E f^P-dpfis) 



4e 



2 ^ Jo dpfdpj 



0<s<t 



jdip^sip^r 

(41) 



/^)-/(p^-£^Ap, 



i=l 



where [p/(.), pj(-)] c denotes the continuous part of [p^(-), /?/(•)]■ 

Let us deal with the case where J ^ 0. If (ej)i<j<p designs the canonical basis of 
R p , then we have p/(t) = (p/,ej) for all t ^ 0. Hence we have dpf{t) = (dpf,ei). As a 
consequence we have for all i G {1, . . . , P} 



P j i(t) = p + f(L{pD,e i )ds+ I ( h k(p J s -),ei)dW k {s) 

+ I I ^k(p J s -), ei}l 0<x<Vk{ p j } [N k (dx, ds) - dxds] . (42) 

, cr Jo JK. 
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It implies that 



[h k (p J s _), et) (h k (p J s ^), a) ds 



E 

fceJUW 

since [Wj(t), = Sijt. Furthermore, if we set by g\{p) = (gk(p),ei), then we get that 

the process 

p 



E 

0<s<t 

-E 

fceJ 



/(rf)-/(p;j-£^^-v/<*) 



i=i 



i=l 



-gi(p J s 



is a martingale. Hence we have 



E 

fceJ 



-E 

fceJ 



i=i 



1 o<x<t, fc (pf„) Ar fc(^,^) 



1 o<x<i, fc (pf_) Ar fc(^ ) ^) 



/(p s j -+^-))-/(p s j -)-E 



-4(p* 



L 0<x<v k (p J s _) 



dxds (43) 



is a martingale because each is a Poisson point process with intensity measure dx ®ds. 
Furthermore we have 

p 



E 

fceJ 

E 



sip 3 ,. + skip',-)) - - J2 ^pU(ftf-) 

i=i ™ 



lo<a;<t) fc (p^_ ) 



fcGJ 



/(p s 7 -+^(p, j -))-/(ps 7 -)-E 



5/(p s j _ 



f [f(pi- + P) - /(P*-) - D p jJ(p)] Tl(p J s _, dp) 



■91(P J S -) 



v k (p J s ^)ds 



(44) 



As the Lebesgue measure of the set of times where p s _ ^ p s is equal to zero, we get that 

f(Pt) - /(Po) - I A J f{p J s _)ds = f{pi)-f{ PQ )- f A J f{ P J s )ds 
Jo Jo 

and it defines a martingale with respect to the natural filtration of (p t ) and the proposition 

is proved. □ 

As announced, the first step consists in proving that equations ( l39i) and ( 140]) admit a 

unique solution. Concerning the equation ( |39l) . we consider the following way of writing 

r t _ p 

P f J = Po 



+ / L(p^)-^^(p g J _)it:eK(p s J _))d S 
P f l f ~ 

+ E / / 9i(Pi-) 1 0<x<Re(v i (^_)) N i{.^d 8 ) i 

i=1 Jo Jr 



(45) 
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and in the same way for equation (HUj) we consider 

fit = Po + f L(p J s _)-J2HpU)Re(v l (pU)ds+ / W^W) 



^-o<x<Re{vi( P J s _))Ni{dx, ds) 

iei Jo Jr 



(46) 



Sufficient conditions (see [18J), in order to prove that equations (l4"5l) and (l4~6l) admit 
a unique solution can be expressed as follows. On the one hand the functions L(.), hi(.) 
and gi(.)Re(vi(.)) must be Lipschitz for all i. On the other hand functions Re(vi(.)) must 
satisfy that there exists a constant K such that we have, for all i and all operator p on 
Ho 

sup \Re(vi(p))\ < K (47) 



Actually such conditions (Lipschitz and (|47p ) are not satisfied by the functions L(.), hi(.), 
Re(vi(.)) and gi(.)Re(vi(.)). However these functions are C°°, hence these conditions are 
in fact locally satisfied. Therefore a truncature method cam be used to make the functions 
L(.), hi{.) and gi(.)Re(vi(.)) Lipschitz and functions Re(vi(.)) bounded. It is described as 
follows. 

Fix k > 0. A truncature method means that we compose the functions L(.), hi(.), 
Re(vi(.) and cji(.)Re(vi(.)) with a truncature function <p k of the form 

4> k {x) = (V>*(ari))i=i,...,p where (48) 
i) k {xi) = -kl Xi <- h + Xil\ Xi \<k + kl Xi ^k (49) 

(50) 

for all x = (xj) G 1R P . Hence, if F is any function defined on M p , we define the function 
F k on R p by 

F\x) = F {<p\x)) 

for all x G M p . By extension we will note F h (p) when we deal with operators on 7i - 

As a consequence, functions L k (.), h k (.) and g k (.)Re(v k (.) become Lipschitz. Further- 
more, as <p k is a bounded function, we have 

sup sup \Re(v k (p))\ < K. 



This theorem follows from these conditions. 

Theorem 3 Let k G M + and let po be any operator on Tio- The following stochastic 
differential equations, in case J = 0, 



Pi 



po + f L k {pi_) - ~g?(pDRe(v k (pU)ds 
J o i= i 

+ Y1 / / aKp'D^xKReivHpi^Niidx^s), (51) 

■ i Jo Jb. 
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and in case J ^ 



Pi = Po+ f L\pl_)-Y,~9-(pi-)R<v^pi_))ds+ fh'M^dW^s) 





(52) 



admit a unique solution. 

Let A k be the infinitesimal generator of the solution of an equation of the form ( [571) or 
. For any f 6 C% and any state p, we have for all k > 1 



where A 3 are the infinitesimal generators defined in Proposition 2. Furthermore in all 
cases, the processes defined by 



where (x) + = max(0,x). 

Proof: The part of this theorem concerning generators is the equivalent of Proposition 
[31 This follows from Proposition [H] and from the fact that, on the set of states S, we have 
4> k (p) = p for all k > 1. Indeed if p = (pi)i=x t ...,p is a state, we have \pi\ < 1 for all i. 

The last part of this theorem follows from properties of Random Poisson Measure 
Theory and is treated in details in [25J for classical Belavkin equations (HI). The proof of 
Theorem 3 follows from Lipschitz character and works of Jacod and Protter in [TH] . 

Let us investigate the proof in the case where J ^ (the case J = is easy to adapt 
to this case with a similar proof). 

Let us prove that equation ( 1521) admits a unique solution (we suppress the index J in 
the solution to lighten the way of writing; we suppress also the index k concerning the 
truncature). As we have sup^ sup peM p \Re(vi(p))\ < K, we can consider Poisson point 
process on R x [0, K]; the part concerning the counting process can be then written as 



Kf{ P )=A 



f{p)- 




(53) 



are counting processes with stochastic intensity 





Hence for all i G J the process 



Ni{t) = card{N i (.,[0,t} x [0,K])} 



(54) 



22 



defines a classical Poisson process of intensity K. As a consequence, for all t, it defines a 
random sequence {(t£, G {1, . . . A/i(£)}} where T l k designs the jump time of A/i(.) and 
the ££'s are independent uniform random variables on [0, -K]. Consequently, the solution 
of the stochastic differential equation is given by 

Pi = Po + L(p s _)ds - ^ / gi(p s -)Re(vi(p s ^)ds 
Jo m Jo 

,t M(t) 

+ E / w^W') + E E ^-'W-f^O (55) 



*eiU{o} iei fc=i 



The solution (1551) is described as follow. Thanks to the Lipschitz property (following from 
the truncature), there exists a unique solution (p\) of the equation 

p\ = p + / L{p\_)ds- E / ft(pi_)-Re(vi(pJ_)ds 
./o ief </o 

+ E / k(pJ-W00 (56) 

For all i E. I, this solution (pj) defines the function i — > Re(vi{p\_)) . We define the random 
stopping time 

Ti = inf 1 t/ E f [ l 0<x<Re ( Vi ( P l_)) N i(dx, ds) > 1 
By definition of Poisson point processes and by independence, we have for all i ^ j: 

P 3t / / lo <x< Re( L ViO^_))Ni(dx,d8)= / / lo<x<Re( Vi (pl_))Nj(dx, ds) =0 

. / Jo J[0,K] Jo J[0,K] 

As a consequence at 7\, there exists a unique index i Tl such that 

i7 M 1,Ml K«J)^ (i ' i!) = 1 ' 

and all the other terms concerning the other Poisson point processes (for different indexes 
of %t x ) are equal to zero. Moreover, we have almost surely 



lo J[o,k] °<*<^K^->J lTlV ' ' ^ "^^KWi-)) 

We define then the solution of floTl) on [0,Ti] in the following way 

(57) 



Pt = p\ on [0,711 



Pn = gi^ipn- 
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The operator p Tl can then be considered as the initial condition of the equation (1551) . 
Therefore we consider for t > T x the process (pf) defined by 




(58) 



ieiU{0} 



In the same fashion as the definition of Ti, we can define the random stopping time T 2 as 



By adapting the expression ( 1571) . we can define the solution on [Ti,T 2 ] and so on. By 
induction, we define then the solution of ( 15T1) . The uniqueness comes from the uniqueness 
of solution for diffusive equations of type of (1551) . Because of the fact that the intensity of 
the counting process is bounded, we do not have time of explosion and we have a solution 
defined for all time t (see (25J or [TSJ for all details concerning such stochastic differential 
equations). □ 

Equations (15T1) or (1521) (with truncature) admit then a unique solution (p t ). In Section 
3, from convergence result, we show that these solutions are valued in the set of states, 
the truncature method will be then not necessary. Therefore solutions of (|51[) and (152"]) 
become solutions of ( 1451) and ( T46l) and as they are valued in the set of states they become 
solutions of ( 1371) and ( 1371) . 

Before to tackle the problem of convergence in Section 3, let us give a proposition 
concerning martingale problem for (A , p ) (f° r some po) an d uniqueness of the solution 
for such problem. This will be namely useful in Section 3. 

Proposition 4 Let po be any operator. Let Af. be the infinitesimal generator of the process 
{p'D, solution of a truncated equation of the form (\51\i or 

The process (p/) is then the unique solution in distribution of the martingale problem 

(A>Po)- 

The fact that the solution of a stochastic differential equation ( 15T1) or ( 1521) is a solution 
of the martingale problem for the corresponding infinitesimal generator follows from Ito 
formula as in Proposition 3. 

This proposition means that all other solution of the martingale problem for (^4 , p ) 
have the same distribution of the solution (p/) of the associated stochastic differential 
equation. This result is classical in Markov Process Generator Theory, it follows from the 
pathwise uniqueness of the solutions of equations ( 15T1) and ( 1521) (see pj] for a complete 
reference about existence and uniqueness of solutions for problems of martingale). 



T 2 = inf \t>T l /J2 I [ 1 



iei Jt ^ J M 



0<x<Re(vi(/? s _))Ni(dx,ds) > 1 
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3 Convergence of Discrete Quantum Trajectories 



In all this section, we consider an observable A of the form ( 1351) with associated subset J 
and / as in Proposition 2. Furthermore we consider an integer k > 1 and the associated 
truncated stochastic differential equations (|5Tj) or (1521) . 

In this section, we show that the discrete quantum trajectory (p J n {t)) (describing the 
successive measurements of A) converges in distribution to the solution of the martingale 

problem for (A k ,po) given by the solution of the corresponding truncated equations (IBTl) 
or (1521) . Next we show that such convergence results allow to conclude that solutions of 
floTT) or fl52l are valued in the set of states. 



Let po be any initial state. In order to prove that the discrete trajectory starting from 
Po converges in distribution, we show at first that the finite dimensional distributions of the 
discrete process converge to the finite dimensional distribution of the solution of the 

martingale problem (A k ,po). Secondly we show that the discrete process (p J n {t)) is tight 
and the convergence follows. For the weak convergence of finite dimensional distributions, 
we use the following theorem of Ethier and Kurtz [TT] translated in the context of quantum 
trajectories. 

Theorem 4 Let A J k be the infinitesimal generator of the solution of the corresponding 
equation (1571) or ( T5B . Let (J 7 ™) be a filtration and let (p J n (.)) be a cadldg adapted- 
process which is relatively compact (or tight). Let po be any state. 
Suppose that: 

1. The martingale problem (Al,po) has a unique solutionfin distribution); 



2. P J n (o) = po; 

3. for all m > 0, for all < U < t 2 < 
and for all f in we have: 



<t m <t<t + s, for all function (9. 



i)i=l, 



lim E 

n— >oo 



f(Pn(t + s))-f(Pn(t) 



t+s 



A J k f{p J n {s))ds 



0. (59) 



Then {p' n {.)) converges in distribution to the solution of the martingale problem for (A J k , po). 

Theorem 4 of Ethier and Kurtz imposes to have uniqueness of solution for the martingale 
problem, this follows from Proposition 4 of Section 2. This Theorem expresses the fact that 
if a subsequence of (p n (t)) converges in distibution to a stochastic process (Y t ), necessarily 
this process is a solution of the problem of martingale associated with (A k ,po).. Indeed, 
from the convergence property (15*9"]) . the process (Y t ) satisfies 



E 



m+s) - f(Y t ) 



t+s 



0. 



(60) 
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As this equality is satisfied for all m > 0, for all < t\ < t 2 < . . . < t m < t < t + s, for all 
function (0j)j=i ... m and for all / in C%, this implies the martingale property of the process 



t - f(Y t ) - f(Y ) 



A J f(Y s )ds. 



Hence, the uniqueness of the solution of the problem of martingale allow to conclude to the 
convergence of finite dimensional distributions and the tightness property allow to conclude 
to the convergence in distribution for stochastic processes. 

Let us deal with the application of Theorem 4 in the context of quantum trajectories. 
Concerning the definition of a filtration (J? 7 ™), we consider the natural filtration of the 
discrete quantum trajectory (p^(t)), that is, if rjn < t < (r + l)/n we have 



J^ = a(p J n (s),s<t) = a(p J p ,p<r). 



It is obvious that JF" 



r /n' 



Assume tightness for instance. In order to conclude, it suffices to prove the assertion 
). As k is supposed to be strictly larger than 1, recall that infinitesimal generators of 
quantum trajectories A J satisfy for all /eC c 2 and for all states p 



The assertion 



A J J(p)=A J f( P ). 
follows then from this proposition. 



Proposition 5 Let po be any state. Let be a quantum trajectory starting from Pq. 

Let Tt be the natural filtration of(p J n (.)). We have: 



lim E 

n^oo 



lim E 

n— too 



f{p J n {t + s))-f{p J n {t) 

f(p J n(t + S)) - f{ P J n (t) 



t+s 



A J k f(p J n (s))ds)l[e t (p^U)) 

' i=l 

\ m 

AJfipi^dsjUeMiu)) 







(61) 



for all m > 0, for all < t\ < ti < . . . < t m < t < t + s, for all functions (0i)j=i,..., m and 
for all f in C%. 

Proof: The discrete quantum trajectory (p J n {t)) is valued in the set of states, we then 
have for all s > 

A J k m{s)) = A j f( P i(s)). 

Let m > 0, let < t\ < t 2 < ■ ■ ■ < t m < t < t + s, let (#1)1=1,...,™ and let / be functions in 
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C%, we have 



E 



E 



E 



f(Pn(t + s))-f( P J n (t) 



t+s 



A J f{pi{s))ds)X[eM(ti)) 



E 



E 



/t+s \ ™- 



f(Pn(t + s))-f(p J n(t) 



t+s 



1=1 



(62) 



Let n be fixed, from definition of infinitesimal generators for Markov chains (see [27]) we 
have that 



fc-i 1 

/fai(*A0) - /(A)) - E -•AZMUM) 



(63) 



i=o 



is a (^T/ n ) martingale (this is the discrete equivalent of solutions for problems of martingale 
for discrete processes). 

Suppose r/n < t < (r + l)/n and i/n < i + s < (Z + l)/n, we have JF™ = J 7 ™^. The 
random states p J n (t) and p^(t + s) satisfy then p J n {t) = p J n {r /n) and p J n (t + s) = p J n {l/n). 
The martingale property ( |63j) implies then 



E 



E\f(pttl/n))-f(pttk/n)/j? /n 



i-i 



E 
E 

+E 



n 

J=k 
■t+s 



A J J(p J n (s))ds/^ 



ft - ^/(^(t)) + (- - (t + .)) ^/(^(t + *)) 



(64) 



As a consequence, we have 



E 



f(p J n (t + s))-f(p J n (t) 



t+s 



t=l 



< E 



^/(A.M)--4 J /(>>;M))A> 



n 

J i=l 



+E 



< Msup 
pes 



[nt] 



n 



Anf( P J n (t)) 



[n(t + s)\ 



A J J(p)-A J f(p) 



L 

H — sup 

n pes 



n 
i J 



(t + s))A J J(p J n (t + s)) 



n 

i=l 
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(65) 



where M and L are constant depending on \\hi\\ and s. Thanks to the condition of uniform 
convergence of Proposition 2, we obtain 



lim 

n— >oo 



E 



Pt+S \ m 

f(p j n(t+s))-f( P j n (t)- / 

^ / <=1 



(66) 



□ 

Finally, in order to apply Theorem 4 of Ethier and Kurtz, it suffices to prove the 
tightness of the discrete quantum trajectory. The following lemma will be useful when we 
deal with this question. 

Lemma 1 There exists a constant Kj such that for all (r, /) e (N*) 2 satisfying r < I, we 
have almost surely: 



E 



< Kj 



n 



(67) 



where = J^ /n = cr{pj,j < r}. 



Proof: Let us deal with the case where J ^ and 1^0 (this is the most general 
case). Let r < I, we have 



E[||pj- P ;|| 2 /A^] =E 



E 



Pl\\ 2 /Mt\ 



Hence we have 



E 



\\ P f-p J r \\ 2 /Mt\ 



E 



E 



E 



j=0 



Y,Vt\pU)-p J r 2 tfipLl) / Mt\ 



.3=0 

p 



E ||4 n) (A-i) - pU + pU - Pr pKpU /m£\ 

.3=0 ' 



4 B) (pf-i) - pU + p/-i - Pi pKpU I M\ n _[ 



(n) 



+ E E 

ie./U{o} 



4 n Vf-i) - pU + pi-i - pi pf(pi-i) / -Mri 



(6£ 



(69) 



As / is supposed to be not empty, for the first term of (16"%]) we have for alH G / 

pKpLi) = l^ipU + o(i)) 
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where functions i>j(.) are defined in Section 2. As C\ n \p) converges uniformly in p 6 S, set 



i? = sup sup sup 

jel n {p,fj,)eS 2 



This constant is finite because all the o's are uniform in p. We have then almost surely 



E E 



4 n) (pf-i) - pf-i + - p, 7 ^ (pf_i) / aci 



< 



card(J) x i? 



n 



For the second term of flUHl) we have 



E E 

E e 

+ E E 

ieJUW 



4 n) (p/-i)-Pz-i P?(p t J -i)/M}-i 



(») 



2i?e ( (4%^) - p^, pf_ x - p r J ) ) pKp/-i) / -M}-i 



<(n) 



+ E E 



\fiLl-fi , r\\'plW-l)/M& 

Concerning the indexes j G </|J{0}' we have 



4 n) (p/-i)-pf-i = ^(^(pf-i)+°(i)) 



n 



In the same way as the constant R, we define 



S= sup supsup \\hj(p) + o(l) || pj(p). 
jeJUW n pes 



For the first term of (170]) . it implies that we have almost surely 



E = 



(n) 



< 



(cardJ +1) x S 



n 



(70) 



Furthermore, as we have ^j 6 jM/ }Pf (pf-i) < 1 almost surely, we get almost surely 



E e 



Itf-l-PrT^-l)/-^! 



< E 



Pf-l-Pr 



-MS 
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For the second term of (1701) . we have 



E E 



2i?e ( (/tV-r) - - p r J > ) p?(pk) / -M^i 



f(«) 



E 



2#e 



( E {^\pU)-pu)m-i),pU-p J M /Mt\ 

WlIM / / 1 



Let us treat this term. As in the proof of Proposition 2 concerning infinitesimal generators, 
with the asymptotic of Cf 1 in this situation, we have uniformly in p G S: 

1 



E (4 l) (p)-p)w'(p) = -(^(p) + °(i)) 



ieJUW 

since the terms in disappear by summing over j £ </|J{0}- As a consequence, by 

defining the finite constant as 



W = sup sup 

n ( PlAt )eS 2 

we have then almost surely 



E b 

j£JU{o} 



2ite[/ E n(^(p)-p)rf 
\ Wll{o} 



(p), P~P 



2Re [(&\pU) - pUpU -pi)) rt{pU) I M™ 



W 
< 

n 



Let us stress that the constant W are independent of I and r. Therefore, we can conclude 
that there exists a constant Kj such that for all r < I, we have almost surely 



E 



\pt-p J X/Mt\ 



<^ + E 

n 



\pU-P j A\ /Mt\ 



It implies that almost surely 



E 



\ P i-p J r \\ 2 /M 



(n) 



<^ + E 

n 

f(n) 



(») 



(71) 



(72) 



Thus by conditioning successively by with i G {2, — r} and by induction, we 

can show 



E 



\pi-P J A 2 /M 



{n) 



< 



Kj{l-r) 



n 



The same results hold when J = or / = by similar computations. 



□ 



This lemma implies the following proposition which expresses the tightness property of 
the discrete quantum trajectory. 
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Proposition 6 Let (p^(£)) be the quantum trajectory. There exists some constant Zj such 
that for all t\ < t < t 2 - 

(73) 



E [\\Pn(t 2 ) - Pn(t)\\ 2 \\ Pn (t) - p n (ti)f] < Zj(t 2 - t,) 2 . 

Therefore, the discrete quantum trajectory {p J n (t)) is tight. 



Proof: The inequality ( 1731) implies the tightness of (p n {t)) (see pa]). Let us prove (1731) . 
It is worth noticing that .M^ = ^T/ n where JF™ is the natural filtration of (p;Q. Thanks 
to the previous lemma we then have: 



< E 



E [||p^(t 2 ) - ^(t)]| 2 ||^(t) - p.(ti)|| 2 ] 
E [IK(KD ~ P J n {[nt\)f\\p J n {[nt\) - P J n {[nt x \)f] 
E [E [\\p J n ([nt 2 ]) - p J n {[nt])f\\p J n {[nt]) - P^KDIlV-^/J] 
E [E [|K(K]) - P J n {[nt])f /Tf nt]/n ] \\ P J n {[nt\) - Pn {[nt x ])f] 
Kj([nt 2 ] - H) 



n 



< 



< 



ATj([nt 2 ] - [nt]) 



|Pn(H) -Pn(K 



E[E [|IPn(K])-Pn(K])ir/^ 1 ]/n]] 



tfj([nt 2 ] - [nt]) #j([nt] - [nh]) 



n 



n 



< Zjfo-h) 2 , 

with Zj = 4(Kj) 2 and the result follows. 



□ 



Therefore we have proved the tightness property of discrete quantum trajectories and 
we can now express the final theorem. 



Theorem 5 Let A be an observable on 71 = C N+1 with spectral decomposition 

v 



ieJIjo 



(74) 



i=0 



where: 



1. For i G {0, . . . ,p} the operators Pi = (Pki)o<k,i<N are the eig en-projectors of A (sat- 
isfying p° 00 ^ 0) 

2. The sets I and J satisfy that I = {i 6 {1, . . . ,p}/p o = 0} an< ^ J = {!> • • • ;P} \ I ■ 

Let p be a state on TCo- Let (p J n (t)) be the discrete quantum trajectory describing the 
repeated quantum measurements of A and starting with po as initial state. 
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1. Suppose J = 0. Then the discrete quantum trajectory {p J n {t)) converges in distribution 
in T>[0,T) for all T to the solution (p J t ) of the stochastic differential equation (\51\i . 
Therefore the process (pf) takes values in the set of states on Tio. The discrete 
quantum trajectory {p J n {t)) converges then in distribution to the unique solution of 
the following jump- diffusion Belavkin equation 

P j t = Po + / L(p J s _)ds + V / / gi(pi-)lo <x<Vlip J_ } [Ni(dx, ds) - dxds] (75) 

J i= i JO Jvi 

2. Suppose J 7^ 0. Then the discrete quantum trajectory (p J n {t)) converges in distribution 
in V[0,T) for all T to the solution (p{) of the stochastic differential equation ff5ij) . 
The process (p/) takes values in the set of states on 7i . The discrete quantum 
trajectory (p J n {t)) converges then in distribution to the unique solution of the following 
jump- diffusion Belavkin equation 

Pi = Po+ I ' L{p J s __)ds+ HpUdWiis) 

^-o<x<Vi(p J _) [Ni(dx, ds) — dxds] . (76) 

ief Jo Jr 
Furthermore the processes defined by 

Nt= [ [ 1 0<x<Mp j ^(dx^s) (77) 
Jo Jr 

are counting processes with stochastic intensities 

/ v i(p J s -)ds. 



Jo 

As in Theorem 3, the last assertion concerning the counting processes of Theorem 5 
follows from properties of Poisson Point processes iVj. It means actually that processes 
defined by 

Vi(p J ,-)ds (78) 



/ / ^x^ipi^NMx.ds) - [ 
Jo Jr Jo 



are martingale with respect to the natural filtration of (pf) (see [],[])• Let us prove the 
convergence results of Theorem 5. 

Proof: In all cases, the convergence result follows from Theorem 4 and Propositions 5 
for the finite dimensional distribution convergence and from proposition 6 for the tightness. 
In order to finish the proof of this theorem, we have to prove that solutions of stochastic 
differential equations ( 15T1) and ( 1521) takes values in the set of states. It is given by the 
convergence in distribution. 
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Indeed, let (p J n {t)) be converging to the corresponding solution (p/) of equation (loTT) or 
we have to prove that this solution is self-adjoint, positive with trace 1. By using the 
convergence in distribution, we have for all zgC 2 

Pn(t) - (p J n(t)Y A pi - {pif (79) 
n— >oo 

Tr[ P m ^ Tr[p-n (80) 

n— *oo 

(z,p J n (t)z) ^ (z,p J t z) (81) 

n—*oo 

where T> denotes the convergence in distribution for processes. As {p J n {t)) takes values in 
the set of states, we have almost surely for all t and all z G C 2 

Pn(*)-(Pn(*))* = 0, Tr\p J t ] = l, (z,p J n (t)z)>0. 

These properties are conserved at the limit in distribution and the process (pj) takes then 
values in the set of states. The proof of Theorem 5 is then complete. □ 



This theorem is then a mathematical and physical justification of stochastic models of 
continuous time quantum measurement theory. Let us stress that in general it is difficult 
to prove that equations () and () admit a unique solution which takes values in the set of 
states. One can notice that such equations preserve the self adjoint and trace properties. 
Concerning the positivity property, it is far from being obvious and it points out the 
importance of the convergence result. 

Let us conclude this article with some remarks concerning these continuous stochastic 
models. 

The first remark concerns the average of solutions of (1751) or (1761) . Let (p{) be a solution 
of f!75|) or fl76l) . In all cases, we have 

E{p{]= f L(E[pi))ds. (82) 
Jo 

This follows namely from martingale property of the Brownian motion and counting pro- 
cesses. This remark concerning (|82|) means that the function 

t - Efan , 

is the solution of the ordinary differential equation 

dp t = L(p t )dt. 

This equation is called the "Master Equation" in quantum mechanics and describes the 
evolution of the reference state of the small system Tio without measurement. In average 
continuous quantum trajectories evolve then as the solution of the Master equation (for all 
measurement experiences) . 
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The second remark concerns the classical Belavkin equations ([3]) and (jl]). In [] and 
[], it is shown that such continuous model are justified from convergence of stochastic 
integral and random coupling method (it does not use infinitesimal generators theory) . 
With Theorem 5, we recover these equations by considering the case where the measured 
observable A is of the form A = XqPq + A\Pi. Indeed in this case, we just have one noise 
at the limit as in the classical case. 

The last remark concerns the uniqueness of a solution of the martingale problems. In 
this article, we have made an identification with the set of operators on T~Lq and R p in order 
to introduce definition of infinitesimal generators and notion of martingale problem (see 
Section 2, Definition 2). As observed, the infinitesimal generators of quantum trajectories 
can be written in term of the partial derivative in the following way 

p 

\ £ Dlf{ hl {p)M P )) = \j2 a M§^ (83) 

p 

D P f{L{p) = Y^h{p) d -M- (84) 



i=l 



by expanding the differential terms. The matrix a(.) = (<%(•) is a semi definite matrix. 
Let W be a P dimensional Brownian motion, the solution of the problem of martingale 
can also be expressed as the solution of 

Pi = Po+ [ L(p J s -)ds+ fa{p J s )dW s 
Jo Jo 

+ Y1 / / 9k(p J s-)l <x<Re(v k (pJ_)) [N k (dx,ds) - dxds] , (85) 
keI Jo Jr 

where a{.) is as matrix defined by cr(.)cr i (.) = a(.) (see []). Let us stress that, in this 
description we deal with a P dimensional Brownian motion corresponding to the dimension 
of R p (which depends only on the dimension of Hq) whereas in Theorem 5 we consider a 
(p + l)-dimensional Brownian motion corresponding to the number of eigenvalues (which 
only depends on the dimension of the interacting quantum system TC). As a consequence 
from uniqueness of martingale problem (Proposition 4) we have two different descriptions 
of continuous quantum trajectories, but they are the same as regards their distributions. 
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